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Abstract 

Application of experimental state data and spectroscopic 
term values shows that the thermodynamic and transport proper- 
ties of reacting PH2-0H2 mixtures are considerably different 
than those of chemically frozen pH2 at temperatures below 300 °R. 
Such differences will be significant from the standpoint of 
nuclear rocket design if it turns out that para- to ortho- 
hydrogen conversion can be catalyzed by radiation in the cold 
dense phase. Calculated H-S data also show that radiation-induced 
PH2-0H2 equilibration at constant enthalpy can produce a tem- 
peiacure drop of at least corresponding to an ideal shaft 

work loss of 15 % or more for a turbine operating downstream from 
the point of conversion. Aside from differences in thermodjmamic 
and transport properties, frozen PH2-0H2 mixtures may differ 
from pure pH2 bn a purely hydrodynamical basis. 


Significance of Possible Para- 
Orthohydrogen Conversion in 
Nuclear Rockets 

Hydrogen is generally a mixture of two 
distinct H2 components. A molecule of the 
component called ”orthohydrogen** may assume 
quantized rotational energies that are not 
accessible to the component called **para- 
hydrogen,** and vice versa. The resulting 
disparities in energy-level populations at 
a given temperature imply differences in 
the thermodynamic and transport properties 
of the two types of hydrogen. Thus, for 
example, the enthalpies and thermal con- 
ductivities of PH2-0H2 mixtures differ ap- 
preciably from those of pure parahydrogen 
in the 36 ^R -500 R range, not only in abso- 
lute magnitude but also in temperature 


dependence. Aside from thermodynamic and 
transport properties, para- and orthohydro- 
gen may differ on a purely hydrodynamical 
basis. The latter possibility, which has 
been given little attention, derives from 
the fact that the differences in accessible 
rotational energies for the two types of 
molecule are associated with differences in 
angular momenta. 

The liquid hydrogen that is stored for 
use in nuclear rocket propulsion is usually 
over 997o pH2* However, the high-tempera- 
ture hydrogen gas that emerges from the noz- 
zle of the nuclear rocket is invariably a 
mixture consisting of 75 % 0H2 and 25 % pH2* 
This change in proportions results from the 
facts that ( 1 ) certain chemical or magnetic 


476 


reactions can convert parahydrogen to 
orthohydrogen when sufficient energy is 
available, and (2) in the presence of a 
suitable conversion catalyst, hydrogen tends 
to a PH2/0H2 ratio that is determined by 
the temperature. Catalysts which convert 
parahydrogen to its equilibrium proportions 
are available in the reactor core (in the 
form of H-atoms from thermal dissociation), 
if not elsewhere. When catalysts are pre- 
sent in large enough quantities that a 
shift in temperature is accompanied, prac- 
tically instantaneously, by a shift in pH2/ 
0H2 ratio to a value corresponding to the 
new temperature, the para-orthohydrogen 
mixture is called '’reacting** or "equilibrat- 
ing.** When no such catalysts are available, 
the PH2/0H2 ratio is not affected by a 
temperature change, and the mixture is 
called **frozen.** In the temperature range 
from.60^R to 300°R the specific heats of 
a given nominal PH2/0H2 mixture differ 
radically according to whether the hydro- 
gen is frozen or reacting. 

In nuclear rocket design, a question a- 
rises as to what type of hydrogen should 
be assumed to exist in the relatively cold 
regions between the propellant tank and 
the core inlet. If the hydrogen is chemi- 
cally frozen throughout this region, then 
it obviously remains in almost pure pHo 
form, and thermodynamic and heat transfer 
calculations can be based on parahydrogen 
properties. If, on the other hand, the 
flowing hydrogen is catalysed by radiation 
or by paramagnetic surfaces, then the hy- 
drogen may be chemically reacting over at 
least part of the region below the reactor 
core, and pure parahydrogen properties may 
not apply. From the standpoint of an over- 
all design study, the effect of catalytic 
conversion is not limited to the values of 
thermodynamic *and transport properties of 
hydrogen at a prescribed temperature and 
pressure. The sudden equilibration of 
hydrogen due to catalysis at a certain 
point in the flow-path can reduce the tem- 
perature by as much as 40^R, through a 
thermochemical process that is analogous 
to the cooling of paramagnetic salts in a 
decreasing magnetic field. Such "conver- 
sion-cooling** can reduce the available 
energy for a turbine operation or alter the 
dimensions of various orifices for an op- 
timum design. 


subject to irradiation and possible en- 
counters with paramagnetic sites, and how 
much conversion can it produce? (2) What 
are the quantitative thermodynamic and 
transport properties of frozen and reacting 
PH2/0H2 mixtures? (3) What is the overall 
effect of conversion on a specific design? 
In regard to the first question, it will 
only be noted that data obtained during 
practical tests indicate that relatively 
’cold parahydrogen can be converted by 
radiation, and that reasonable chemical 
models have been advanced to explain this 
effect. The third question is in the pro- 
vince of rocket designers. The results 
given here pertain only to the second ques- 
tion, viz., information on the basic pro- 
perties of PH2-0H2 mixtures. 

Equilibrium Orthohydrogen 
Fraction 


The equilibrium orthohydrogen fraction 
at a given temperature is independent of 
pressure, at least for pressures in the 
range of nuclear rocket application and 
lower. For temperatures below 500^ (900%) 
the equilibritjm orthohydrogen fraction, 
is given to an accuracy of 0.01% by 


^ (odd) .0 /x/® -0./T 

We “ ^ ^ 8i ® 

j-1 ^ / j-0 


( 1 ) 


where j is the number of a rotational level. 
Even and odd values of j correspond to pH2 
and 0H2 rotational levels, respectively. 

The parameter 9x represents the molecular 
rotational energy of a particular level 
divided by Boltzmann’s constant. Using the 
spectroscopic term values cited in the 
definitive work of Woolley, et al. at NBS 
(Ref. 1), it is found that a suitable ap- 
proximation is 

0j(®K) - 85.3757 j(j+l) - 0.065769 

+ 0.0000712 (2) 


for values of j less than or equal to 8. 
This prescription ignores the role of vi- 
brational excitation, which is negligible 
below SOO^K. The parameter gi in Equation 
(1) is a statistical weight given by 

gj - 3(2j+l), j odd 


Whether hydrogen conversion in the sub- 
core region is of design significance de- 
pends on three questions. These are: (1) 

does an effective mechanism for conversion 
actually exist in relatively cold hydrogen 


- (2j+l), j even (3) 

At temperatures above T ■ 500°K, the equi- 
librium orthohydrogen fraction is ■ 
0.750. 
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Evaluation of Equation (1) yields the de- 
pendence of equilibrium orthohydrogen frac- 
tion on temperature shown in Figure 1. It 
is seen that the temperature range over 
which cJq varies significantly is 30°K to 
220°K, or 54°R to 396°R. 

Calculation of Enthalpy and Entropy 

The calculation of enthalpies and entro- 
pies is obviated by expressing the results 

as ideal gas values, augmented by "real- 
gas** terms that arise from intermolecular 
reactions. The latter are obtained from 
an equation of state, which is usually 
written in the form of a temperature and 
density dependent for the compressibility 
factor Z ■ pv/RT. In obtaining the pre- 
sent results it has been assumed that the 
equation of state 

Z - Z(T,P) (4) 

is the same for either para- or ortho- 
hydrogen. The state equation used is that 
given by Wooley et al. for normal hydrogen 
(75% ortho). This state equation, which 
was based on rather old normal hydrogen 
data, has been compared with data reported 
more recently by Michels et al. (Ref. 2). 
T)rpical comparisons are shown in Figure 2. 
An analysis of the difference between the 
old and new data, which is taken as a 
measure of the present uncertainty in the 
state equation for normal hydrogen, shows 
that the corresponding uncertainties in 
enthalpy and energy of normal hydrogen are 
respectively only about 1 Btu/lb and 0.002 
Btu/lb °R. Furthermore, recent data by 
Hermans et al. (Ref. 3), on the difference 
between the second virial coefficients of 
para- and normal hydrogen, indicate that 
the use of a normal hydrogen state equa- 
tion on PH 2 - 0 H 2 mixtures will not lead to 
an error greater than that corresponding 
to the present uncertainty in the normal- 
hydrogen state data. 

The real-gas enthalpies and entropies 
of a number of mixtures have been calcu- 
lated using tabular data given by Woolley 
et al. to obtain the integrals involving Z. 
Ideal-gas values of the enthalpy and en- 
tropy for both frozen and reacting PH 2 - 0 H 2 
mixtures have been obtained from the 0^ and 
gj defined in Equations (2) and (3), ^ 
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The enthalpies of frozen pH2, frozen 0 H 2 
and equilibrating hydrogen at 0 and 800 
psi are shown in Figure 3. The correspond- 
ing specific heats, obtained by differen- 
tiation of H with respect to temperature, 
are shown in Figure 4. It is seen that 
EH 2 behaves as a mixture with an interme- 
diate specific heat only above 100°K; at 
lower temperatures the specific heat of 
reacting hydrogen is appreciably higher 
than that of either frozen component. By 
contrast, the specific heat of a frozen 
mixture is obtained by simply calculating 
the average of the frozen pH 2 and 0 H 2 val- 
ues, weighted by their respective concen- 
trations in the mixture. 


478 






The viscosity of a PH 2 - 0 H 2 mixture is 
determined by transport of momentum alone, 
and is therefore independent of the pH2- 
0 H 2 composition insofar as pH 2 and 0 H 2 con- 
form to the same equation of state. On 
the other hand, the thermal conductivity of 
an ideal-gas mixture of para- and ortho- 
hydrogen is related to the viscosity, 77 , 
by Eucken's relation 



Kq - (97-5)0° 7?/(4M) (5) 

where Cv is the specific heat of the ideal 
gas, y is the ratio of specific heats, 
and M is the molecular weight. The real- 
gas thermal conductivity K is related to 
Kq by the following Enskog relation: 

K/Kq = 1 + 0.575E + 0.5017E^ - 0.204e3 (6) 

where E is given by 

E = Z - 1 + T (||)p (7) 

Figure 5 shows the thermal conductivity 
for pHo, 0 H 2 and a reacting mixture at 0 
psi and 800 psi. As in the case of speci- 
fic heat, the thermal conductivity depends 
strongly on whether a given mixture is 
frozen or reacting. 


Cooling Due to Conversion 


If a chemically frozen parahydrogen gas 
is suddenly exposed to a conversion-catalyst 
such as radiation or a paramagnetic sur- 
face, the resulting equilibration to a 
PH 2 - 0 H 2 mixture is an irreversible process. 
Hence, the reaction cannot be isentropic 
even if it happens to be adiabatic. In 
general, the enthalpy change due to a dif- 
ferential conversion process is 


dH = TdS + Vdp + li^^dx^ + /i2^^2 

where x^^ and X 2 are the pH 2 and 0 H 2 mole 
fractions, respectively, and where and 
/i 2 are the chemical potentials. In the 
case where no heat enters or leaves the 
system and the pressure is kept at a con- 
stant value. Equation (8) simplifies to 


TdS = (/7^-/i2)dx2 


(9) 
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This increase in entropy is accompanied by 
a decrease in temperature, which can, for 
example, subsequently reduce the amount of 
available energy from a turbine operating 
between two prescribed pressures. The 
possible effect of such cooling is seen 
from the approximate H-S chart shown in 
Figure 6, which has been constructed from 

Figure 3 and calculated entropy data. 

When parahydrogen at 297°R is converted to 
607o 0 H 2 at constant enthalpy and constant 
pressure, the temperature decreases from 
297^R to 269^. If conversion had not 
occurred, the energy available to a tur- 
bine operating between 975 psi and 675 psi 
would be 111 Btu/lb. However, it is seen 
that 60% conversion would reduce the avail- 
able energy to 98 Btu/lb. The magnitudes 
of such losses depend upon the inlet en- 
thalpy and, of course, the amount of con- 
version that is assumed. Figure 7 shows 
the percent loss in ideal energy for an 
assumed conversion of 607® as a function of 
inlet enthalpy. 

Possible Hydrodynamic Differences 
Between pH p and oH^ 

The difference between the thermodynamic 
and transport properties of parahydrogen 
and those of orthohydrogen are completely 
specified by the rotational energy levels 
accessible to each, at least insofar as 
the two modifications conform to the same 
equation of state. However, from the hy- 
drodynamical standpoint, another type of 
difference might be worth exploring. 

Parahydrogen molecules may take intrin- 
sic angular momenta of 0, 2h, 4h, etc.; 
whereas orthohydrogen molecules may take 
intrinsic angular momenta of h, 3h, 5h, 
etc. At temperatures below 180 k, over 
977o of all parahydrogen molecules are in 
a state of zero angular momentum. For 
practical purposes one may say that, at 
these temperatures, parahydrogen is a 
fluid composed of molecules with no intrin- 
sic angular momentum, whereas orthohydro- 
gen is a fluid whose molecules have an 
intrinsic angular momentum of fi. Hence, 
on a rigorous hydrod 3 mamical basis (Ref. 4), 
the equations for energy and momentum con- 
servation, which are sufficient to describe 
the behavior of a parahydrogen fluid, must 
be augmented by an equation of angular 
momentum in the case of orthohydrogen. In 
the language of transport theory, colli- 
sions between ortbohydrogen molecules in- 
volve an extra summational invariant. It 
is therefore conceivable that the details 


of such phenomena as the transition from 
laminar to turbulent flow might differ in 
the two types of fluids, independently of 
any effects associated with thermodynamic 
and transport properties. 
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